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BULK DIFFUSION IN A SYSTEM WITH SITE DISORDER 

By Jeremy Quastel 1 

University of Toronto 

We consider a system of random walks in a random environment 
interacting via exclusion. The model is reversible with respect to a 
family of disordered Bernoulli measures. Assuming some weak mixing 
conditions, it is shown that, under diffusive scaling, the system has a 
deterministic hydrodynamic limit which holds for almost every real- 
ization of the environment. The limit is a nonlinear diffusion equation 
with diffusion coefficient given by a variational formula. The model 
is nongradient and the method used is the "long jump" variation of 
the standard nongradient method, which is a type of renormalization. 
The proof is valid in all dimensions. 

1. Introduction. Consider a system of particles occupying sites of a mul- 
tidimensional integer lattice. The particles are attempting jumps to nearest 
neighbor sites at rates which depend on both their position and the objec- 
tive site. The rates themselves come from a quenched random field, and are 
chosen so that the system satisfies a detailed balance condition with respect 
to a family of random Bernoulli measures (the random field Ising model 
at infinite temperature). The interaction between the particles is given by 
a hard core exclusion rule: Attempted jumps to occupied sites are simply 
suppressed. 

Such systems have been used to model electron transport in doped crystals 
[1, 2, 3, 8, 10, 12, 14]. In this case, the hard core exclusion rule is given by 
the Pauli exclusion principle. The crystal itself creates a periodic field in 
which the particles move. In the presence of impurities, the field is random. 

The purpose of this article is to study the transport properties of such 
a system. In particular we are interested in the influence of the random 
field on the rate of bulk diffusion. In the absence of the exclusion rule, 
this is simply the result of diffusive scaling of a single particle moving with 
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reversible rates in the random field. At the other extreme, if the exclusion 
is present but the field is constant, the bulk diffusion also turns out to be 
independent of the density of particles. However, when both the random 
field and the exclusion rule are present, one sees a nonlinear dependence of 
the bulk diffusion on the density. This diffusion coefficient can be computed 
by the Green-Kubo formula, but we have not been able to find such a 
computation in the literature. We give a variational formula for the bulk 
diffusion, which is equivalent to the Green-Kubo formula, and which we can 
establish rigorously. 

We should point out that in principle it is not even clear that under 
diffusive scaling such a system has a hydrodynamic limit. The main work 
here is to show that in fact all the influences of the random field are, on a 
sufficiently large scale, contained in a diffusion coefficient, which depends on 
the statistics of the field, but not on the randomness itself. One of the main 
steps in this direction is to show that the system in a box of side length L 
has a spectral gap no smaller than cL~ 2 . 

To establish the hydrodynamic limit, one needs to prove some version 
of Fick's law, namely, to replace the microscopic current by the gradient 
of the density field multiplied by the diffusion coefficient. The system we 
are considering turns out to be of nongradient type. Roughly speaking, the 
gradient condition says that the microscopic current is already of gradient 
form. At this time, the only method for the nongradient systems is the one 
developed by Varadhan. The idea is to replace the current by a gradient plus 
a fluctuation term. However, in our case, such a decomposition cannot hold 
microscopically, because the fluctuations of the gradient of the density field 
arising from the random field are very large, and one can only make sense of 
these gradients over sufficiently large mesoscopic distances where stochastic 
fluctuations are reduced by the central limit theorem. It is necessary to 
perform the decomposition in some mesoscopic scale. 

The model we are considering has been studied in the physics literature 
by means of rough approximations and Monte Carlo simulations. The case 
of a period two field in one dimension has been solved exactly [20, 24]. 
Surprisingly, the resulting equation is linear. This special case turns out 
to be an "almost gradient system," by which we mean that the fluctuation 
term is explicit. In other words, one can find the minimizer in the variational 
problem for the diffusion coefficient. This is certainly not the case for other 
fields, and it can be shown that the diffusion coefficient is nonconstant in 
general. 

The main result of this paper was described in detail in [16] along with 
a sketch of the proof and an unpublished, incomplete manuscript [18]. An 
important motivation for the method described there is that the nongradi- 
ent method as developed in [17, 22, 23] does not work in low dimensions. In 
higher dimensions, such an approach can be made to work by subtracting 



BULK DIFFUSION IN A SYSTEM WITH SITE DISORDER 



3 



a term from the microscopic gradient term to make it mean zero. One then 
has the nontrivial problem to show that the subtracted term vanishes in the 
limit, and this can be done if d > 3. This was suggested in [16], and a proof 
following these lines appeared recently [6]. This inspired us to write up the 
details of our unpublished manuscript carefully as the present paper. The 
advantage of the original approach in [16] is that it works in all dimensions 
and elucidates the connection to renormalization. The present article is in- 
dependent of [6]; the only result we have used beyond [16] and [18] is the 
proof of the spectral gap of the Bernoulli-Laplace version of the model by 
[4], which improves on what was written in [18]. Note that both [4] and [6] 
need the moving particle lemma (Lemmas 5.1 and 5.2 of this article), which 
is in some sense what tells you that the system is diffusive. 

In addition, we prove here the continuity of the diffusion coefficient in the 
full interval [0, 1] and the uniqueness of the hydrodynamic equation. The 
uniqueness needed to be added as an assumption in [6], and the continuity 
was proven only on (0, 1). Note that the continuity of the diffusion coefficient 
at is of considerable interest because there it becomes the asymptotic 
diffusivity of a single particle in the medium, a classical homogenization 
problem with a different variational formula that has to be related to ours. 

The review article [16] contains a sketch of the main idea of the proof 
contained in this article. The sketch contains a small misprint which has 
unfortunately led to a great deal of confusion. The formula for a certain 
central limit theorem variance, 

(1.1) limsupE[(i^~ 1 Wx, {-K 2 L K y l K~ l W K )\ = Cm 2 {l - m 2 )<7~\ 

K^oo 

is stated without a numerical factor C on the right-hand side. We would 
like to explain the nature of this minor mistake, which does not affect the 
integrity of the proof. In fact, as long as it is finite, the exact value of 
the constant C is not even relevant in the proof. To explain this, we need 
to describe the basic idea of the proof. This is a bit easier in one space 
dimension. 

The main problem in proving the hydrodynamic limit is to understand the 
large scale behavior of the average current w XtX +i between site x and x + 1, 
multiplied by a factor e , the size of the system. The current is a function 
of the particle numbers and field values at those two sites. One is allowed 
to take some average of these currents over a larger, but still microscopic 
block. Write the sum on the block as w Xy y — ^/z= x ^z,z+i- Averaging does 
not appear to help too much unless w x ,x+i were of gradient form, r x +ih— r x h 
for some h. Then there would be an easy term by term cancellation. But 
in this model the current is not of gradient form unless a = 0. It might 
work, however, if we could find an h such that the difference can be made 
appropriately small. This is the standard approach to nongradient systems, 
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which as explained, works for this model only in higher dimensions. Our 
alternative approach is to use a long jump version of the current, w Xj y, which 
is the same function evaluated with variables at sites x and y, but now x and 
y are chosen far apart. It turns out to be true that one can replace w x ,y by 
Cw Xt y, where C is related to the diffusion coefficient. This helps somewhat 
because one needs to normalize the sum in w x ,y, and a normalization factor 
of y — x in the denominator is now available. But to compensate for the 
huge e _1 factor in front of the current, one would have to take y — x = e _1 
which is no longer microscopic. So one needs something more. One can get 
the needed extra factor by averaging the long jump ends x and y over some 
large microscopic blocks, say, x E Ai and y E A2, where Ai and A2 are two 
blocks of side length K beside each other. Then one gets a free factor K~ 1 / 2 
{K~ d l 2 in dimension d) from the central limit theorem, and together with 
the K -1 from the typical distance between x and y, this gives a prefactor 
j£-(d+2)/2 w \ l [ c \ l can heat the e~ l if K is chosen as large as e~ 2 /( d+2 \ 
So we want to replace 

2/-1 

(1.2) K^Av^Y, w 2 ,2+i by if 1 Wk — K 1 Av x& A 1 w Xj y. 

y£A 2 z=x y£A 2 

If one instead tries to replace Av ZtZ+ i e \ lU \ 2 w Z)Z+ i by K~ 1 Wk = K -1 x 

^WreAi w X)V , the resulting central limit theorem variance will not vanish. 

S/SA 2 

This is not a surprise and is simply because one forgot the extra average 
on the left-hand side of (1.2), so the corresponding variables are weighted 
differently than those on the right-hand side of (1.2). But it appears some 
readers drew the erroneous conclusion that this is because of the missing 
factor C in (1.1). We regret the misprint and the confusion it has caused, 
but it does not disqualify in any way the main idea of the proof as described 
in [16] and [18]. 

The proof of the hydrodynamic limit contained in this article is essentially 
the same as that of [16] and [18]. The main new material is the proof of the 
continuity of the diffusion coefficient, the uniqueness of the hydrodynamic 
equation and the proof of (1.1), the details of which were not included in 
either [16] or [18]. This is in Sections 7 and 8. It requires a nontrivial ex- 
tension of the standard nongradient computations, as one is missing in this 
problem the usual average over translations. 

The paper is organized as follows. In Section 2 we describe the model 
and the main results. In Section 3 we give a sketch of the proof of the hy- 
drodynamic limit, emphasizing the new problems that arise because of the 
random field. Sections 4-8 contain the details of the proof. In Section 4 we 
recall some standard results in perturbation theory which allow us to use 
the variance method for nongradient systems on functions whose range is 
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up to a small constant times e~ 2 /( rf+2 ) as long as they satisfy an "integra- 
tion by parts" estimate. In Section 4 we prove the key "moving particles 
lemma" which gives the spectral gap and the two-block estimate and the 
"integration by parts" estimate for the long jump current. Sections 6-10 are 
devoted to the computation of asymptotic variances which is the heart of 
the nongradient method. The standard material is in Sections 6 and 7. In 
Sections 8 and 9 the computation is extended to the "long jump current" 
where one no longer has an average over shifts. In Section 10 we complete 
the proof by showing that the long jump current can be replaced by its 
average with respect to local equilibrium at the scale e -2 /( rf + 2 ) where the 
fluctuations have been sufficiently dampened. Section 11 contains the results 
about continuity and a type of Holder continuity of the diffusion coefficient 
and Section 12 the uniqueness of the hydrodynamic equation. 

2. The model. Let a x , x G Z d , be a bounded (\a x \ < B) stationary, 
ergodic random field satisfying the following mixing conditions: For some 
7 > max{4, 2(d + 2)/d 2 }, there is a constant C < oo such that, for all I > 0, 
for all / G T Al with E[f] = 0, 

(2.1) E[\Av xeAK r x ff) < c(K-Hr d ' 2 nri 

where t x is the shift by x, Ak is the box of side length K, and E is ex- 
pectation with respect to the field. Typical examples are a periodic, or a x 
could be independent, identically distributed random variables. But (2.1) is 
quite general. For example, one could divide Z, d into boxes of side length A. 
On each of the boxes one could pick at random from a finite list of patterns 
(functions on {1, . . . ,A} d ). 

For each e = L -1 , L a positive integer, we have a system of N = 0(e~ d ) 
particles on eL d j7L d moving in this field. At most one particle is allowed at 
each site. A particle at x attempts to jump to nearest neighbor sites y at 
rate 

(2.2) e~ 2 (l + e a y- a *). 

If there is no particle in the way, the particle is allowed to jump. However, if 
there is a particle in the way, the jump is suppressed, and everything starts 
again. All the particles are doing this independently of each other, and since 
time is continuous, one can ignore the occasion of two particles trying to 
jump onto each other simultaneously. 

The state space of our process is {0, l} eZ / z . Configurations are denoted 
r\. r] x = 1 or depending on whether there is or is not a particle at x. Our 
system is a Markov process on this state space with a generator given by 
e~ 2 L £ , where 

(2.3) L £ = ^ a xy(v)^xy, 

\x-y\=l 
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where the rates are given by 

(2.4) a xy (r t ) = l + e- ( - a y- a * )( > r, y-Vx) ) 

and the lattice gradient is given by 

(2-5) V xy f(v) = f(T xy v)-f(v), 

where T xy r\ represents r/ with the occupation numbers at x and y exchanged. 
For each A, the generator is reversible (self-adjoint) with respect to the 
product measure 

(2.6) Z _1 exp ^2 (ot x + X)r) x . 



.re 



,1 



Z is the normalization. The parameter A, which is called the chemical po- 
tential, can be adjusted to vary the average density of particles. The relation 
between A and the density, m, is as follows: For < m < 1, A is chosen so 
that 



(2.7) E 



1 + e a +X 



m. 



We could alternatively fix the number of particles so that the density was 
m. Then the process is reversible and ergodic with respect to that measure 

(2.8) Z^ 1 cxp^ 

a xVx 3{Av x ri x =rn} ■ 

Taken over the allowable m = ie d , i = 0, . . . , e~ d , these give us a full set of 
ergodic invariant measures. Note that the measures (2.8) are not simple to 
describe. The Dirichlet form is given by 

(2.9) D(f) = E[f, (-L)' l f) = £ E[{V xy ff). 

x~y 

The empirical density field [i e € M(T d ), the set of measures on the d- 
dimensional torus T d , is given by 

(2.10) ^ £ {d9) = e d Y J Vx 'ex j 

x 

where 5q gives mass one to the point 6 T d . Let us choose initial dis- 
tributions of our process so that, for some given smooth m° :T d — > [0,1], 
/i £ =^ m° (9) d6 in probability. The corresponding probability measure on 
Z)([0,T] — > M(T d )), the space of right continuous trajectories with left limits 
in M(T d ), will be denoted P £ . Of course, P £ depend on the field, a. 
We can now state our main result. 
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Theorem 1. For almost every realization a of the random field, P e 



*m(t,6)d8> as e 

(2.11) 



dm 

— — = V • D(m)Vm, 
dt y ' 



0, where m is the unique weak solution of 

m(O,0) = m°(0). 

The diffusion matrix D{m) is a nonrandom continuous function on [0, 1]. It 
is given by the following formulae. For m G (0, 1), 

D(m) = o~(m)A' (m), 

where the conductivity, a, is a symmetric matrix whose associated quadratic 
form is given by the variational formula, 



(2.12) (PMm)P) 



2infE 

g 



\e=l 



The infimum is taken over all local functions g(r], a) of the configuration 
and the field. The shift is given by T x g(rj,a) = g(r x r],T x a). The expectation 
E is over the random field a and the expectation (-) m is over the infinite 
product measure (2.6), where A is chosen as in (2.7). For m = 0, D(0) is 
the limiting covariance of a free particle in the field, given by the classical 
homogenization formula, 



(2.13) (J3,D(O)0) 



2z~ x inf E 



J2(e ao + e ae )(P e + T~ e U - Uf 



where z = E[e Q °]. For m = 1, D(l) is the limiting covariance of a test "hole" 
which coincides with the limiting covariance of a free particle in the field 
a = —a. 



The result was obtained earlier in [6] in d > 3 and under the assumption 
that the diffusion coefficient is continuous and that the limiting equation 
has a unique solution. 

A few comments follow. 

Regularity of D(m). From (2.12), a{m) is upper semicontinuous for m G 
(0, 1). It is not hard to check that A'(m) is continuous on (0, 1) and, therefore, 
D{m) is upper semicontinuous on (0,1). The test function g = in (2.12) 
shows that a{m) < Cm(l — m) for some C < oo and it is elementary to 
check that \'{m) < C/m(l — m) for another finite C. Hence, D(m) < cl for 
some c < oo for all m £ (0, 1). From the moving particles Lemma 5.2, one 
obtains also a lower bound D(m) > c -1 / for m G (0, 1). In Section 9 we show 
that D{m) is Holder 1/2 with a coefficient which may behave badly at the 
edge; \D(m\) — D(m2)\ 2 < C(mi(l — mi))~ l \mi —m2\ for some C < oo. This 
follows from the characterization of the diffusion coefficient in Sections 6 and 
7 without too much work. It is also shown there that D(m) is continuous 
on the whole interval [0, 1]. 
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Weak solutions. By a weak solution of (2.11), we mean a function m : [0, T] x 



T — ► [0, 1] satisfying 



(2.14) 



iVrn' 2 
o Jjd m( 



l —d8dt = suv\ / [mV- <j>- \<p\ 2 m(l -m)]dOdt 

1 — m) <$, Uo Jf d 



< oo 



for each T < oo and for smooth test functions 6 on T d x [0,T], 
t=T 



(2.15) 



4>m d6 



. T [ ^-mdOdt- [ T [ Vd>- D(m)Vmd6dt. 

*=0 Jo JT d Ot Jo JT d 



In Section 10 we show that, under the Holder continuity proved in Section 9, 
such weak solutions are unique. 



Associated dynamics. One can produce other dynamics with (2.8) as er- 
godic reversible measures. Consider the dynamics associated to the Dirichlet 
form 



(2.16) 



where p is finite range and symmetric and b X;X+y (r]) = T x bo ty (r]) for some 
finite range &o,y( ? ?) bounded above and below (for y in the range of p). The 
expectation is with respect to any one of the measures (2.8). The corre- 
sponding dynamics have a particle at x attempting to jump to x + y at rate 
Py(boy( T xV) + bo y (To y T x r])e ax+y ~ ax ). Since our methods are based on esti- 
mates involving the Dirichlet form, they extend easily to equivalent Dirich- 
let forms, and one obtains an identical theorem with (2.12) replaced by the 
infimum over 



(2.17) 



2E 



X!m b °y \P ' y^y - m) - v o,s/ J2 T *9 



A natural example is &oe = e~ ao Vo(^ — Ve) + e _ae ?? e (l — t/o). The resulting 
dynamics have a particle at x attempting to jump to each nearest neighbor 
site at rate e~ ax . 



Mixing conditions. The mixing conditions given in (2.1) are chosen for 
convenience and are not meant to be optimal. They are mainly to point 
out that the theorem holds under some very weak mixing conditions on the 
variables a x . A nice problem is to consider the case of unbounded field a 
(the condition that \a x \ < B is very important for the method). 
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3. Hydrodynamic limit. In this section we give a sketch of the proof 
of the hydrodynamic limit. Many of the arguments are now standard and 
can be found, for example, in [11]. However, at some points new ideas are 
needed, especially in low dimensions where there is not sufficient averaging to 
control the fluctuations from the random field. We will sketch the approach, 
emphasizing where new methods are needed, leaving the rigorous proofs for 
Sections 4-8. 

The evolution of the empirical density field \x e is described by the following 
set of stochastic integral equations: 

(■T 



4>(t,e)fi £ (t,d8) 



t=T 



t=0 



T d dt 



(t,0)n e (t,dO)dt 



/V 2 £ m,ey) 
Jo , . 



<f>(t,£x))w xy dt 



\x-y\=l 



+ f T e d £ m,ey) 
Jo ii. 



4>(t,ex))V xy ridM, 



xy 



\x-y\=l 



Here M xy are independent Poisson "sawtooth" martingales running at rates 

Ixy 

(3.1) V) a 



')-%(l-%)(l + e a ^). 

J' E[Lf 2 - 2fLf) ds if M f (T) = f(r,(T))- 
/(t?(0))- ft Lf( V (s))ds, true for any Markov process, one easily computes 
the quadratic variation of the martingale term, 

2i 



Jxy — wxy{Vy ~ Vx) — — VxjO- + e 

Using the formula E[\M f (T)\ 21 
rT 



(3.2) E 1 



I e d ((f>(t,£y)-(f>(t,ex))V xy ridM xy 
\x-y\=l i 



<Ce a 



where C depends only on T and £ d J2xeeZ d /z d <p 2 (x). 

If one starts with a nondegenerate product invariant measure (-E[??o] 7^ 
or 1), we have the equilibrium, or stationary process which we denote by Q e 
and from the bound on the entropy of any initial distribution with respect 
to that reference measure, we obtain directly the bounds 

(3.3) H(P £ /Q £ ) < Ce~ d , f D(^J t )dt<Ce 2 ~ d , 

Jo 

where ft is the marginal density of the nonequilibrium process at time t. If 
V is any bounded function, we can estimate by the Feynman-Kac formula, 
the spectral theorem for self-adjoint operators, and the variational formula 
for the principle eigenvalue of L + V, 



T" 1 \osE Qf - 



cxp 



V(r)(s))ds 
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(3.4) 



< sn J(Vf)-e- 2 £ ((V^v//) 2 )} 

I \x-y\=l > 



\x-y\-- 

where the expectations are with respect to the invariant measure and the 
supremum is over densities relative to that measure. The superexponential 
bound 



(3.5) 



expe d ^ 4>(ex) f w xx+e (s)ds 
<2exp(t-s) Y 4> 2 ( x ) 



x£el 

is obtained this way using the integration by parts formula 
(3-6) (w X yf) =-^({Vy-Vx)V xy f), 

true for any / and any invariant measure on any set containing x and y. By 
Schwarz's inequality, we have 

(3-7) (w x , x+e f)<C((V xy ^f)Y 2 , 

which gives (3.5). Once one has (3.2) and (3.5), tightness of the measures P e 
follows by a standard argument from Garsia's lemma (see [11] for details). 

Once one has tightness, it remains to identify the limit measure and for 
this, we take the limit of the stochastic integral equation. From (3.2), the 
martingale term is asymptotically trivial. Hence, the work is to identify the 
limit of the term involving the current in terms of the empirical density field. 
Since the sites are distance e apart, 

e _1 (^(ey) - <j>(ex)) = (y - x) ■ V<j)(ex) + o(l). 

Fix a direction eo and call J = eo • V0. Our job is to identify the limit of 

(3.8) f T e d J{ex)e- l w x , x+eo dt 

for smooth functions J on the T d as 

(3.9) [ T [ J(9)y2D eoe (m(9,t))d e m(9,t)dt, 
Jo Jr d g 

where m(6,t) is the density of the limit of the empirical density fi e . It is 
relatively easy to see that the latter is well approximated by something of 
the form 
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where fj x = Av y&A ir\y and A x is a cube of side length £ around x, since 

m(t,9)d9 is the weak limit of fj, e (t,d9) and fj s x £ (t) = fj, £ (t, {9 : \0 - x\ <5}). 

In order to make this replacement, one could try to replace (3.8) by some- 
thing like 

f T 

(3.11) / e d ^2 J{£x)e~ l T x Avy & \ K w y ^ + e dt 1 



where Ak is a box of side length K about the origin. Note that it is not 
hard to make such a replacement. Performing a summation by parts on the 
difference between (3.8) and (3.11), one obtains an error of J^Tdt, where 

(3.12) T = e d Yl [J(ex)-Av yeAK J(e(x + y))}e" 1 



One easily estimates from (3.7) 

(3.13) (17) < CK e ||VJ|| oo ^ eZd/£ - lzd£ - 1 ((V x , a;+e0 y7) 2 ) 1/2 , 
so that 

(3.14) (17) - £ d ~ 2 D(^f) < C'K 2 e 2 

for some new C < oo. So as long as K = o(e _1 ), such a replacement can be 
performed. On the other hand, it is not so clear how replacing the current 
by its average really helps us to get closer to something like (3.10) or where 
the nontrivial term a(m) would come from. 

Instead we will choose some functions Q e K which depend on variables in 
box Ak of side length K, large with e _1 , so that the latter is well approxi- 
mated by 

(3.15) ( T e d J2 J{ex)e- l T x ® e K dt. 

x& d je-^TA,e 

The functions Q K are given explicitly [see (3.28)], but for now we do not 
need the explicit form to explain the basic argument. 

There are also a class of objects for which one can readily check the 
asymptotics are trivial. Let g(r], a) be any local function and define the shift 
T x g = g(r x r], T x a). Then 
r T 

(3.16) / e d J2 J{ex)e~ l LT x gdt 

x& d /e- 1 1 d 

is asymptotically trivial. This can be seen by Ito's formula, which says that 
the above is equal to 

rT 



£ d+i J(ex)T x g T - [ e d+1 V J(ex)a X)X+e V X)X+e g 

Jo 



dM xx+e . 
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Terms such as (3.16) are called fluctuation terms. Let 

(3-17) Vk = Av xeA > K w x , x+eo - Lr x g - ^ lsr> 

e 

where A' K is a box of side length K — £, where I is chosen so that Vk depends 
only on variables in A' K . We need to show that 



(3.18) 



/ e d J(ex)£~ 1 t x Vk ds ^ in P £ probability. 

From either the entropy bound in (3.3) together with the entropy inequality, 
or directly from the Dirichlet form bound in (3.3) to show (3.18) for a 
function Vk, it suffices to prove 

(3.19) inflimsupsup{ e ^ e " 1 ( T ^/)- £ " 2 E(( V ^+eV / 7) 2 )! , = 0- 

9 f[x e J 

Taking the sum out of the supremum, we obtain an upper bound of the form 
inf limsupsup< e~ 1 {Av x& x K VkJ) ~ e~ 2 Av x& A K Ei^+^j 2 ) (■ 

9 e-fO [ e J 

The expectation is now over a canonical measure (fixed density of particles) 
on the box of side length K and the supremum is over all relative 
density functions, as well as all densities. Letting Lk denote the generator 
corresponding to the Dirichlet form J2xeA K Ee((^+eV7) 2 )i we recognize 
this as the variational formula for 

where Xk, e is the principle eigenvalue of 

K 2 L K + eK d+2 V K . 

Note that K 2 Lk is used because it has a spectral gap of order one. If we 
write down the formal Rayleigh-Schrodinger series for Ak,e, we find 

£ -2 K - {d+ 2) XK£ 

= s- 2 K^ d+2 \eK d+2 (V K ) + e 2 K 2( - d+2 \V K , (—K 2 Lk)~ 1 Vk) + ■•■}• 

The functions Q K are specially chosen so that they have mean zero with 
respect to any canonical measure on A^-- This is also true of the currents 
and the fluctuation terms and, hence, (Vk) = 0. Therefore, 



(3.20) f e d J2 J{ex)r x [e~ l V K -a K ]ds 
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is asymptotically trivial where 

(3.21) a K = K d+2 (V K ,(-K 2 L K )- 1 V K ). 

In Sections 4 and 5 this is proved for K < coe~ 2 ^ d+2 \ where Co is a small 
constant as long as we have an estimate of the form 

(3.22) {V K f) 2 <CK d J2 <(V*,v7) 2 > 

\x-y\=l,x,y£A K 

holding for some C < oo for all densities / on Ak- For our special choice of 
Off [see (3.28)], (3.22) will follow from (3.7) and the moving particles lemma 
proved in Section 5. We want K as large as possible to control fluctuations 
from the random field and, hence, we will always choose 

(3.23) K = c e- 2 ^ d+2 \ 
Let 

(3.24) a K = a K -E[a K ]. 

We claim that, for any fixed g, almost surely in the random field, 

(3.25) lime d ^ J(ex)T x a K = 

in P e probability. To see this, note that ax, which is a function of m = re- 
taking values in iK~ d , i = 0, . . . , K d , can be easily extended to a continuous 
function of m G [0, 1] by linear interpolation. For each 5 with <5 _1 a positive 
integer, divide Z rf /e _1 Z d into disjoint boxes of side length fe _1 and label 
them (3. Let f)p be the particle density in box (3. By the two block estimate, 

(3.26) lim \ime d Av Av x£ pJ(ex)[a K {f}x,K,r x a.) - a K (fjp,T x a)] = 

o — »0 s — >U p 

in P £ probability. Now ax (jn) is a function of a x , x € Ak, hence, by the 
mixing condition (2.1), for 7 > 1, 



E 



Av Av xe pJ(ex)a K (r)p,T x a) 



71 

< A-vE[\Av xe f3J(ex)a K (Vi3,T x a)\' y } 



= 0([5eKy d / 2 ). 

For fixed 5, this is 0{e ld ' V 2 ( d+2 )) which is summable in e^ 1 = 1,2,... as 
long as 7 > 2(d + 2)/d 2 . By Chebyshev's inequality and the Borel-Cantelli 
lemma, the term goes to zero for almost every realization of the random 
field. This proves (3.25). 

Hence, we have reduced the problem to proving that 

(3.27) inf lim supE[ax (m,g)] = 

9 K— >oo m 
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for our specific choice of &x, as well as proving that (3.10) is well approxi- 
mated by (3.15). 

The choice of functions Q e K is not unique. The standard choice is some- 
thing of the form DV e n, where V e r/ = rj e — i]o, making the passage from (3.15) 
to (3.10) easy. However, V e ?y does not have mean zero, so one has to sub- 
tract a term E[V e r]\f]i] for some large I, and try to deal with that term in 
a different way. One can check the size of the subtracted term after appro- 
priate averaging is only small in dimensions three or higher. In [17] it was 
suggested that the standard approach could work in d > 3 and it was carried 
out in [6]. 

We choose instead in (3.15) 

(3.28) Q e K = "(VA K )Wft, 
where, for any integer I, 

(3.29) Wf = r 1 w At)TuA t 

is the block renormalized or long jump current. Here A £ = {1, . . . ,£} d and 
for any two nonintersecting subsets A and B of Z d , 

(3.30) W A ,B = Av X £ A ,y£BW xy 

is the average current over A and B, where w xy are given by (3.1). The 
prefactor v is given by 

(3.31) v(m) = cj eoe (m) /m(l — m). 

Computations in Section 8 explain why v has to have this form. (3.27) is 
proved in Sections 6 and 7. 

4. Perturbation theory. We recall some standard results from perturba- 
tion theory which we will be using in a specific context. 

Lemma 4.1. Let Hq be a nonnegative self-adjoint operator on a Hilbert 
space with Ao = a simple eigenvalue, and spectral gap Ai > 1. Let V be a 
real potential bounded by 

(4.1) (u,Vu) <a((u,H u) + l/4(u,u)) 

for some a < 2/5. Let H = Hq + V . Then for |A| = 1/2, the Green function 
(A — H)~ l has the convergent series 

oo 

(4.2) (A - Hy 1 = £ [(A - Hoy'vnX - H y\ 

n=0 

with 

(4.3) ||[(A - Ho^VriX - Hor'W < 5(5a/2)". 
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PROOF. Note that (H + 1/4)~ 1/2 is well defined. By (4.1), 
(4.4) (u, (H + l/4)- 1 / 2 F(F + l/^)~ l/2 u) < a\\u\\l 

Since V is real, (H + l/4y 1/2 V(H + 1/4) -1 / 2 i s self-adjoint, and therefore 

(4.4) implies the operator bound 

(4.5) \\(H + l/±r 1/2 V(H + 1/4)- 1 / 2 !! < a. 
Rewrite [(A - Ho)- l V] n as 

(H + \y 1/2 [(Ho + \)(X - Ho)-\Ho + \T 1/2 V(Ho + \y 1/2 ] n 

(4.6) 

x(#o + |) 1/2 . 

Since |A| = 1/2 and Hq has gap greater than one, we have ||(i?o + l/4) _1 ^ 2 || < 
2, \\(H + 1/4)(A - J^o)" 1 !! < 5/2, ||(flo + 1/4) 1 /2( A - tf )-i|| < 5/2. This 
proves the bound and the convergence. □ 

The spectral projection for the ground state of H is given by 

(4.7) P = ^-l ^-H)- 1 - 

2m J\\\=i/2 

Hence, we have a convergent expansion for P. Using this expression, one 
obtains the familiar Rayleigh-Schrodinger series (see [19]) for the ground 
state energy of H + V with the nth term bounded by 5(5a/2) n . 

Corollary 4.2. Let W be a real potential and 

(4.8) H = -K 2 L K + eK d+2 W. 

If K 2 Lk has gap of order one, then the Rayleigh-Schrodinger series for the 
ground state energy converges provided W satisfies 

(4.9) (u,Wu)<K- d ^D K (u) 1/2 \\u\\ 2 , K <c e- 2,{d+2) 
for some Co small enough, or 

(4.10) Halloo <Ci, K< (C7ie/10)- 1/(d+2) . 
Furthermore, in both cases one has 

(4.11) e~ 2 K ~ d ~ 2 infspecF < K d (W, -(L K )~ l W) + o(l). 

Proof. Assume the gap of K 2 Lk is one. From (4.9) and Schwarz's 
inequality, 



(4.12) eK d+2 (u,Wu) <eK^ d+2 ^ 2 [K 2 D K (u) + l/A\\ 



un|] 
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By Lemma 4.1, the power series for the Green function converges. Further- 
more, the nth term is bounded by c(eK^ d+2 ^ 2 ) n . Therefore, the nth term 
in the Rayleigh-Schrodinger series is bounded by c(eK^ d+2 ^ 2 ) n . One can 
compute that the first term is zero from the assumption (4.9) on W. The 
second term gives 

(4.13) {e 2 K d+2 )K d (W,-(L K )- 1 W}, 
while the other terms are bounded by 

(4.14) ( e id d + 2 )/ 2 ) 3 . 

Now suppose that (4.10) holds instead. We can choose a in the previous 
lemma to be eK d+2 Ci. The first two terms in the Rayleigh-Schrodinger 
series can be computed as before. The other terms are again bounded by 
o(l) after multiplying e~ 2 K~ d ~ 2 . This proves the corollary. □ 

5. Moving particle lemma. We need a preliminary result. 

Lemma 5.1. Suppose p is a homogeneous Bernoulli measure (a x = 0) 
on Z 1 conditioned to have N particles. Let k be a positive integer and p x a 
sequence of positive numbers with Ylx=i = 1- Then 

k—l 

(5-1) E,[(V hk f) 2 ] <J2p^E,[(V x , x+ if) 2 ]. 

x=l 

It is important that there is no multiplicative constant on the right-hand 
side. 



Proof of Lemma 5.1. We prove it by induction. For k = 3, it is ele- 
mentary to check the inequality directly. Suppose it is true for k — l and let 
p\ H h Pk = 1- Let q = p2 + • ■ • + Pk ■ From the k = 3 case we have 

(5.2) £[(V 1)fe /) 2 ] < p^EKV^f) 2 ] + q- l E[(V^ k f) 2 ]. 

From the inductive hypothesis, since q~ 1 p2 H + q~ l pk = L 

fe-l 

(5.3) ^[(V 2>fc /) 3 ]<gX;p^[(V^ +1 /) 2 ]. 

x=2 

This completes the proof. □ 

Lemma 5.2 (Moving particle lemma). Suppose p is a Bernoulli measure 
on Z 1 with external field a taking values in [—B,+B]. Then 

(5.4) [\f(T 1L n)-m] 2 drtr,)<e 13B L £ [ (V x , x+1 f) 2 dp,( V ). 

J l<x<L-l J 
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Proof. Suppose we change the measure \x to a new measure jl by chang- 
ing each on to the nearest value of the form Bj/L, j an integer. Since there 
is always such a point with |a« — Bj/L\ < B(2L)~ 1 , the Radon-Nikodym 
derivative dfi/d/l is bounded above and below uniformly by e B l 2 and e~ B l 2 , 
respectively. Therefore, at the cost of a factor of e B , we may assume that a 
takes values in {Bj/L : j = —L, . . . , L}. By the same reasoning, at the price 
of a factor e 4jB , we may assume that ao = ai = B. 

Let A be the set 

A = {xi-.a Xi =K,i = l,...,k}. 

By definition, 

1 l,L r ) — J-x^,X2 ' ' ' 1 x k _ 2 ,x k _ 1 1 x k ,x k _ 1 ' " ' I x 3 ,x 2 1 x 2 ,x 1 r l- 

By Lemma 5.1 with p" 1 



X s + l-X s ' 



WPU*?) - /(t?)) 2 ] < e 4 ^ E ~ ^T ^[(/(r W ii7) " /(»?)) 

s= l x s+l x s 



2i 



We have to bound 



E,[(f(Tx s ,x a+1 v)-f(v)f] 



x s+ i - X 

We are now in the same situation as before except no a x can take value K 
when x s < x < x s+ %. Let us change a Xs and a Xa+1 to the value B(L — 1)/L. 
The price we pay is a factor exp{25L -1 }. Continuing this procedure, we 
have a proof of the lemma. □ 



From the moving particles lemma, we obtain as a consequence two basic 
bounds, the spectral gap and the two block estimate. The arguments leading 
from Lemma 5.2 to the two block estimate are completely standard (see 
[11]), so we simply state the result we need. We denote by fj x the empirical 
density of particles on a box of side length n around x. Also {-) m denotes 
the expectation with respect to the fixed total density of particles fj E -i = m. 



Lemma 5.3 (Two block estimate). Let < K < oo and a any field with 
—K <a x < K . Let F be continuous function on [0, 1] x [0, 1] and J a smooth 
test function on T d . Let K — > oo as e — > with K<5e~ 1 . Let 

(5.5) F £ , j = e ^JM(F(4,4)-^r l ,r 1 )). 

X 
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Then for every realization of the field a with —K < a x < K , 



lim sup lim sup sup sup< (y e> sf)j 

8->0 e->°o \y\,\z\<K f { 
0<m<l 



E ((v^ +e y/) 2 ) m <o. 



(5.6) 



The expectation (-) m is with respect to the canonical measure (2.8) with fixed 
density m on Z d /e~ lr L d and the supremum is over all relative densities f . 

The following result from [4] gives the spectral gap of the Bernoulli- 
Laplace version of our model to the correct order. 

Lemma 5.4. Let < K < oo. There exists a constant C = C{K) < oo 
such that, for any field a with —K < a x < K , any A, and any < N < \ A\, 
for any f : {0, 1} A M, 



(5.7) 



Var A) jv(/) < CIA]" 1 ^ 



v 



.x,y£A 

Here VarA,Ar,_ft" and E\^ t x denote the variance and expectation with respect 
to the measure in our random field on A conditioned to have N particles. 

Together with the moving particles lemma, one obtains the spectral gap 
of the nearest neighbor dynamics to the correct order. 

Theorem 2 (Spectral gap). For each K > 0, there exists a C < oo such 
that, for all a with —K < a x < K , all cubes A^ of side length L, all < 
N < L d , and all f : {0, 1} Al -» R, 



(5.8) 



Var Aii7V (/) < CL 2 E Al 



Proof. By Lemma 5.4, we have 



Var Ai ,7v(/) < 



C 

A7 



Ea l .n 



E (v,,/) 2 

x,y£k L 
\x-y\=l 



E ( v ^/) s 

x,y&K L 



For each x,y E A^, choose a canonical path x = x\,X2, ■ ■ ■ ,x n = y with Xi 
and Xj+i by moving first in the first coordinate direction, then in the second 
coordinate direction, and so on. By the moving particles lemma, we have 



EA L , N [U(TxyV)-f(v)Y]<e 



13K 



n 



E E ALjN [(V Xi , Xi+ jf]. 



Ki<n-1 
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Summing over x and y, noting that n < dL and that each nearest neighbor 
pair is used for the path between d(L/2) d+1 pairs x and y, we obtain the 
result. □ 




6. The diffusion coefficient. We now discuss the computation of the 
asymptotic variance (3.27). Suppose f{rj,a) is a function depending on the 
particle and field configuration in some finite box and which has mean 
zero with respect to every invariant measure for our process on that box. 
We denote the class of such functions Ge- Denote by Q the increasing union 
of Ge- For any v G Ge, we form the sum 

(6.1) V K = Y, t * v > 

x£A' k 

where A' K is a box of side length K — £ and we define 

(6.2) V K = (V K , (-Lk^Vk) = sup j 2(V K f) - 



The expectation is with respect to an extremal invariant measure, that is, 
the canonical invariant measure (2.8) with fixed number of particles, and the 
supremum is over all densities with respect to such a measure. In particular, 
\ K depends on fj\ K = Av x ^\ K ri x , as well as a x , x G Ak, 

(6.3) \ K = V K (VA K ,ai). 

Since v G Ge, if K ^> £, we can write v = L^h for some h supported on 
Ak- Then Vk = J2xeA' K T x^K^ = J2xeA' K T xh an d, hence, we have V k < 
CK d . For each fixed m G (0, 1), define 

(6.4) [v,v] =\imsupK~ d E[V K (m K ,a)]. 

Although [v,v] is well defined for each v EG, it is not so clear how to 
compute it. For this purpose, we introduce an auxiliary Hilbert space TC and 
compute [v,v] for objects in G by mapping them to TL. Let 7r m denote the 
product measure (2.6) with density m. A form £ e is a function of r] and a 
depending on the basic basis elements e of Z, d . We make a Hilbert space Ti 
out of them through the inner product 



(6-5) « e ,o> = £E[<e 



e)m]- 



For g a local function, consider the exact form = J2y£Z d ^bT y g indexed 
over nearest neighbor bonds b = x,x + e in Z d . Since g is local, there are 
only finitely many nonzero terms so the sum is well defined. Note that 
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is covariant in the sense that T x £b = £, Tx b- So they can all be reconstructed 
out of the basic forms £ e , e running over basis elements of This defines 
a Hilbert subspace E of exact covariant forms. An exact form possesses the 
algebraic property of closedness: If b\, bi, ■ ■ ■ , b n are an ordered set of bonds 
making a loop in Z d , then J27=iCb i ( r ] bl bi ~ 1 ) = 0. An example of a form 
which is closed but not exact is given by (Vr/) e = V e r/ = n e — t]q. Let C be 
the closure in 7i of the closed forms. One can check by standard arguments 
(see [11]) that E has codimension d in C and that C is the closure of the 
linear span of E and V77. 

We now compute [v, v] for some special cases. When v = Lg for some local 
g, we use the map 

(6.6) g i-> £ e = V , e ^2 T *9 

x 

and find 

(6-7) [Lg,Lg] = {(£,£)). 

When v = iuo,e , we have the image w 1— > Vrj and 

(6.8) [w,w] = ((Vri,Vri)). 
Furthermore, for v = Lg, the map above gives 

(6.9) Ke,u] = -|E[(V^,e e >m]. 
The variational formula (2.12) for a translates to 



(6.10) /3-<tP = w£ 



PeWQ,e ~ Lg 



Fix a direction e. For each fixed £, Wf given by (3.29) is embedded in £, as 
wf. From (3.7) and the moving particles lemma, one has a uniform bound 
on the norm [W/,W/]. Furthermore, one can check by a straightforward 
computation that, for any g £ Q, 

(6.11) }im[Wf,Lg}=0 
and less trivially that 

(6.12) lim [(3 -W t ,(3- W t ] = m 2 (l - mf (3 ■ a' 1 p. 

I— too 

In other words, the limit W of the Wt represents an (unnormalized) orthog- 
onal complement to Lg with respect to [•,■]. From the above computations, 
we have 

(6.13) lim \ni\{w- Lg-uWA 2 } =0, 

1^00 9 
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where v(m) = a(m) /m(l — m), which is the first step in the renormalization. 

To prove (6.12), we take the limit in variational formula (6.2) (see [11]) 
for any v G Qi to find 



(6.14) 



[v,v] = sup 



1 2 ]T E[(V e T_ x V,T_ x ^ 



>]-«£,£» , 



where —L\ e V = v and £ is a closed form. Since £ is closed, it can be approx- 
imated by /3Vry + V J2y T y9- The first term can be computed explicitly and 
the variational formula becomes 



(6.15) 



[v,v] = sup { 2E 
9e<5,/3 1 



;^,-vi> 9 ) 2 ))}. 



Applying this to Wg gives (6.12). It is worth noting that 



E 



and 



(3Vr] 




[v,(3-w-Lg] 



[03 -w- Lgf 



so this also shows that Q is a Hilbert space and Q = w © LQ. 

At this point it is worth making a quick remark about associated dy- 
namics. If one chooses instead a dynamics as in (2.16), then one will have 
a corresponding generator L and a corresponding current w different from 
the w above, but satisfying all the needed estimates. One checks in ex- 
actly the same way as above that lim^oo inf 3 [(w — Lg — DWf) 2 ] = 0, where 
P(m) = a(m)/m(l — m) and a is given by (2.17). The key point is that the 
choice Wi is independent of the particular choice of model. This is important 
in later sections where the special form of Wg, is used repeatedly. 



7. Structure of the gradient space. The state space on which we work 
is ft = [—B,B] Z x {0, 1} Z points of which we call (a,rj). Fix < m < 1. 
On {0, l} zd , we have the product measure 

(7.1) n = Z~ 1 expl ^ (a x + A(m))% >, 



22 



J. QUASTEL 



with A(m) chosen as in (2.7). Denote corresponding expectation by .£>«[•]■ 
On £1 we have the corresponding annealed measure 

(7.2) E[F(a, V )]=E[E a [F]). 

A collection of functions {u>b} where b runs over unoriented bonds of 7L d 
is called a translation covariant form if, for all bonds b and sites x in Z d , 

(7.3) T x uj b = uj Txb . 

Here the shift operator t x acts on functions F by (r x F)(a,rj) = F(T x a,T x r]), 
where (r x v)y = Vx+y and {T x a) y = ct x + y . A translation covariant form is rep- 
resented by {w e }, where e runs over the basic bonds e\ = (1, 0, . . . ,0), . . . , = 
(0, . . . , 0, 1) since every b = r x e for some x and e and ujb = T x uj e . The collec- 
tion of translation covariant forms is then a Hilbert space Ti with the norm 
||K}|| = £ e E[ We 2 ]. 

An ordered finite set of bonds b\ , . . . , b n is called a closed loop if T bri ■ ■ ■ T bl r] = 
7} for all 7]. A form {cob} is called closed if for any closed loop of bonds 
bi,...,b n , 

n 

(7-4) J2 UJ h(a,T bi _ 1 ---T h r ] ) = 0. 

i=l 

A translation covariant form is called exact if oj e = J2 x ^ei~ x g for some local 
function g(a,rj). Note that only finitely many terms in the sum are finite, 
so this makes sense. Note also that an exact form is automatically closed. 

Denote by C the closure in Ti of the space of square integrable translation 
covariant closed forms and by £ the closure of those that are exact. We have 
£ C C and, in fact, the containment is strict as V e r/, the form (r/(i 5 o,...,o) — 
^(o.—.o)) ■ • • > V(o,. ..,o,i) — ^(o,...^)) i s closed but not in the closure of the span of 
exact forms. The main result follows: 

Theorem 3. C = Vrj®£. 

Note that the theorem has been proved in different situations in [6, 17, 22]. 
The proof here is very similar. It requires that the disorder field a be a 
stationary ergodic process, but does not use the mixing conditions (2.1). 

PROOF of Theorem 3. Let u>b be a closed form. Let Ak be a box in 
Z d of side length K centered at the origin. For any set A, let r]\,a\ denote 
the collection of variables rj x ,a x , x S A. For bonds b inside Ak, let 

(7.5) = ~E[u b \a AK ,r]A K ]- 

ujf satisfies the closedness condition (7.4) for closed loops bi,...,b n con- 
tained inside Ak - The set {0, 1} Ak is divided into ergodic classes J2xeA K Vx = 
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iV and on each we can sum along bonds to produce an unambiguous function 
Sk with 

(7.6) V b S K = u? 

for any bond b inside A#- . There is a free variable depending on J2xga k Vx 
and, therefore, we can choose Sk so that E^ISkI^xgAkVx] = as well. 
Let 

(7.7) G K = E[S 3K \aA K ,r] AK ] 
and, for any bond b in Z rf , 

(7.8) tf = K~ d J2V b T x G K . 

Note that £^ makes sense because only finitely many terms are nonzero, 
and that it is by definition an exact form. If e is a basic bond, we can write 
£g = £l ,K + , where £,\' K corresponds to terms in the sum with x and 
x + e in Kk and £ 2,K corresponds to terms in the sum with one of x and 
x + e in Ax and one not in Ak ■ 

il' K = K~ d T- x V X:X+e E[S 3K \a y ,r]y,yeA K } 

x,x+e£An 

(7-9) =K~ d E[r_ x V x , s+ e5 3 A'|Q!r_.A JC ,^_ B A J r] 

x,x+e£Apc 

= K~ d H^e\a T ^ x A K ,Vr. x A K ]- 

x,x+e£Apc 

By the martingale convergence theorem, we have — > w e as K — > oo. 
Note that 

(7-10) #* = cfV e »7, 

where does not depend on 770 or r/ e . Suppose we take a bond 6 which 
does not have any vertex in common with e: 

(7.11) v b & K = K- d J2 V e r x V Txb G K 

\{x,x+e}C\A K \=l 

SO 

(7.12) E[|V 6 ^| 2 ]<CK-( d+1 ) n\V x+b G K \ 2 ]<C'K- 2 . 

\{x,x+e}nA K \=l 

In other words, E[(V b cf ) 2 ] -> as if — ► 00. Since Cg does not depend on 
7/0 or r/ e and this is true for any b not having a vertex in common with e, 
any limit c e of the c^f cannot depend on the configuration 77. 
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We have thus shown that if u>b is a closed form, then there exist exact 
functions ^ and a function such that ^ — — ► u> e and any weak 
limit £g of the £g'^ must be of the form 

(7.13) c e (a)V e r ? . 

Assume for a moment that such a weak limit exists. We now want to show 
that a form of the type (7.13) is in V e f? © £ as well. 

First of all, we can always subtract cV e i] from (7.13), so taking c = 
E[c e (a)], we can assume without loss of generality that E[c e ] = 0. 

Now we solve explicitly for a function Gk such that V&Gr- = c e (T x cn)^/ bV 
for bonds b inside Ak of the form (x,x + e). Trying Gk = Y^xeA K a x{ot)r] x , 
one obtains the set of equations 

(7.14) a x - a y = c e (T x a). 

These are readily solved by taking ao = and for any other x £ Ak , 

n 

(7.15) a x = -J2c ei {T Xi a), 

i=l 

where (xi,Xi + ei), i = 1, .. . , n, is a sequence of bonds from to x. It does 
not depend on the path taken because (7.13) is closed. 

Let £g be as in (7.8) with this new Gk, and again break it up as £,f = 
£,l ,K + fe' • The first term Q K converges to (7.13) in the same way as 
before. The difference is that we now compute S^' K directly. It is given by 

(7.16) g'* = rfV e »7, r« = K~ d £ T- X a x 

plus an analogous term with x + e G Ax and x ^ Ax- Connect boundary 
points x of Ak to the origin in some deterministic way: Say, (xj(x),Xj(x) + 
ei(x)), i = 1, . . . , n(x), with n(x) < ifd. Then 

n(x) 

(7-17) r? = -K~ d ]T E c ei(s)( T ^(s)-* a )- 

i£A,3:+e^Ajf i=l 

As -RT — > oo, this converges to E[c e ] = by the ergodic theorem. This proves 
that forms which can be written as in (7.13) are in V e r](B£. 

Finally, we have to prove the key analytic point which is the boundedness 
in norm of the defined after (7.8) and, hence, the existence of a weak 
limit point. This is in fact the hard point in the nongradient method and 
the crucial point where the spectral gap is used. We will use the following 
lemma. 
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Lemma 7.1. Let Ai and A2 be disjoint finite subsets ofL d and x a point 
not in either. There exists a constant C < 00 such that, for all functions f 
ofr]AiUA 2 u{x}> all a and al1 V G A 2; 



(7.18) 



E a [(V xy E a [f\r] AlU{x} ]) 2 \r] Al ] 

KCi^^Ealif-EMMft+Av^EaliVxzf) 2 ]^]}. 



Proof. The proof is standard. We write E[-] for -Eo[ - |?7Ai] and P for 
the corresponding probabilities. We can assume without loss of generality 
that E[f]=0. Note that 

(7.19) (V xy E[f\ V{x} }) 2 = (E[f\i lx = 1] - E[f\ Vx = 0}ft( Vx jt Vy ) 
and 

E[f\r ]x = l]-E[f\r ]x = 0] 

= E[f\ Vx = 1] - E[f\r) x = 1,^ = 0]+ Elfin* = hVz = 0] 
- E[f\r] x = 0, Vz = l}+ E[f\ Vx = 0, Vz = l]- E[f\ Vx = 0] 



(7.20) 



E 



/I 



1 - ffe 



Pz 



Tfe = 1 



+ E 



1 - nz 



Pz 



nx = 1 



/1 



1 -ft 



??x = 



where p z = -Ef^]- Since < m < 1 and \a x \ < B < 00, we have a bound 
8 < P2 < 1 — <5 for some <5 > 0. By Schwarz's inequality, we therefore have 

(E[f\r ]x = l]-E[f\r ]x = 0]) 2 

(7.21) 

< C{E[f 2 \i lx = 1] + £[(V^/) 2 |7fe = 1] + £[/ 2 |? fe = 0]}. 

But again E[\q x = 1] < C£[-] and E[-|t7 x = 0] < OE[-] and the lemma fol- 
lows. □ 

Applying the lemma with / = E a [S3K\VA K ], x G A^, y = x + e ^ A#-, 
Ai = Aj<- \ {x}, A2 = {x + e, x + 2e, . . . , x + Ke}, we obtain 

(7.22) E a [(V x , x+e f) 2 ] < C{K-^ a [S 3 ^; Sax] + Av zeA2 E a [(V xz S 3K ) 2 }}. 

By the moving particles lemma, Lemma 5.2, 

z 

E a [(y xzSzk) 2 ] <CK^^E a [{y x+ ^i) ejX+ie S 3 K) ] 



(7.23) 



i=l 



n(z) 



CKj2E a [(u : 



3K v 
x+(i—l)e,x+ie' 



i=l 
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By the spectral gap, Theorem 2, there is a C < oo such that, for any a, 
(7.24) E a [S 2 K ]<CK 2 

beA K 



Taking expectation of (7.22), (7.23) and (7.24) over a, and using Jensen's 
inequality and E [u;&] = C < oo , 

(7.25) E[(V,, x+e G^) 2 ] < CK 2 

and, hence, another application of Jensen's inequality gives the required 
bound on the norm of £, 2,K : 



(7.26) E[\g K \ 2 ]=E 



2n 



\ one of 2,2i+e£Ajf / 



< c. 

□ 



8. The long jump current. So far everything we have described is stan- 
dard (see [11]). Applying the perturbation theory and the discussion from 
Section 3, we have obtained that the difference between (3.8) and 

(8.1) e d J{ex)T x e- l v e ^ e (r] kK )Av y ^ K T y Wtdt 

goes to in P e -probability as e — ► followed by I — > oo, where K = cq£~ 2 ^ ( q!+2 ) . 
However, we will need to take the £ in w\ on a much larger scale in order for 
the averaging to beat the fluctuations from the random field. We will again 
use central limit variance computation as in the previous section, but now 
we will not be allowed to have the full average over shifts which is crucial in 

(6.2) , so the problem has to be handled in a new way. 

To compute the central limit theorem variance of the long jump current, 
we will use its precise form and a renormalization procedure. 
Using its precise form, we can rewrite the current (3.1) as 

(8.2) Wxy =riy-rjx + £,xCy ~ (x£,y, 
where 

(8.3) & = (1-Vx)e a ', (x = Vxe~ a *. 

Now we introduce some notation. Let Q be a subset of Z d containing two 
nonintersecting subsets Ai and A2. The average current on Ai,A2 is given 
by 

WA U A 2 = Av X £A lt y<zA 2 W xy . 

Let 



77A = AvzcAVx 
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be the empirical density on a set A C Z d and A a = AA(m) be the empirical 
chemical potential defined implicitly through 

(8.4) m = Av x( za{i] x ); = Av x£ a » • 

We will often write Aa for Aa(^a)- Define also 

(8.5) £a = Av xeA £x, Ca = Av x€A (x- 
We have 

(8-6) w Ai ,a 2 = % 2 - %x + CaiCa 2 - Cai6v 2 - 

The corrected average current is given by 

(8-7) W Al ,A 2 (. m ) = W Al ,A 2 - [7A 2 -7Ai], 

where 

7A = 2m(l - m)A A + me'^iU - (Ca)~ Xa ) 



(l-m)e x ^HU-(U)x A )- 



Note that wa 1 a 2 is a function which depends on the variables a x and r\ x 
for x G Ai U A2, as well as an additional variable m G [0, 1]. The form of the 
correction is chosen so that there is a C < 00 such that if Ai and A2 are 
cubes of side length £, for all m G [0, 1], all ficZ^ containing Ai U A2, 

(8.9) E[E[{w Al)A2 {m)} 2 \m = m}} < Cr 2d . 

We will prove this in the next section as Lemma 9.1. 
We now set up the renormalization procedure. 

Fix a large positive integer M and let be the union of the cube 
A\j = {1, . . . , M} d and its translate A 2 M = TMeA l M . 
Let rj, i = l,..., M, be defined through 

(8-10) ±Av£ l Av]*? M+1 J2 at = Av™rm 

k=i 

for all dj, i = 1,...,M. Although it is elementary to write down a closed 
form expression for the rj, we will only ever use (8.10). Next define 

(8.11) p x =r x . e , x£fl M . 

Note that Av xe Q M p x = 1 and the average Av x ^n M p x a x of a x , x G £Im has 
the special property that if a x = A x+e — A x , then 

(8.12) Av xe a M p x a x = —[Av yeA 2 M A y - Av xeA i M A x ]. 
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Lemma 8.1. Let M = K/£ be an integer. Recall the definition (3.29) 
long jump current 

wf=r 1 w A t tTuA t. 

Let 

(8.13) V K/ = Av x ^ K/lPx T lx W!-W e K . 

There is a C < oo such that, for all m £ [0, 1], and positive integers M and 

e, 

(8.14) # d E[£[^(-L Q J~V^fe K = m]} < CM d+2 £~ d . 

Proof. From (8.12), 

£~ 1 Av x( zn K/e p x Ti x [T ee jA - 7a] 

(8.15) 

= K 1 [Av y€A 2 K/ T ey j A - Av x€A i^T ex j A ] , 
so we can write Av ye n K/e p x T£ x Wg - Wfc = A K>i + B K>i , where 

(8.16) A K ,i = Av xe o, K/l p x n x i~ 1 w Aijn<iAv 

(8.17) B K>e = Av x&k i K/ ^ y&A 2 K/ K~ l w TlxAlyTtyAl . 

The w are defined in (8.7). By Jensen's inequality, 

(8.18) E[E[A 2 K/ \fjn K = m}) < Av x£ n K/e p x r 2 T ex E[E[w 2 \fjn K ]] < C£~ 2 ~ 2d . 
In the same way, 

(8.19) E[E[B 2 K/ \fj UK = m}] < CK~ 2 r 2d . 
By spectral gap, 

(8.20) E[V K {-Ln K r l V K \m K ] < CK 2 E[Vl\fjn K }. 

We conclude that the left-hand side of (8.14) is bounded by CK d £~ 2d . Since 
K > £, this is again bounded by CK d+2 £~ 2d ~ 2 , which is the same as the 
right-hand side of (8.14). □ 

Now let K = £M N for some N and £ n = £M n . R# n+1 is readily seen to 
be a union \J xG q i~£ nX Qi n of copies of Qg n . So a point x £ &£ n+1 can be 
represented as (xi,^), where x\ G Qi n and X2 G £Im- Continuing in this 
way, we have 

and ffl d nQ J{ -nl 
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If x £ fi]{f, we can write x = (x±, . . . ,x n ), where X\ € Q-m and if a x is a 
function on fi^, we can define the average 

(8.21) Av xeQ n i a x p'£ l > = Av Xie n M ■ ■ ■ Av Xn£ n M a ixu _ tXn) p Xl ■■■p Xn . 
The main result of this section is the following: 

Theorem 4. Fix M a large integer and K = £M N . Let 

(8.22) R Kjt = Av xeadnUK r x WipW - w e K . 
Then 

(8.23) lim Jim K d E[E[R K/ , {-Ln K R K ,e\fjn K }} = 0. 

I— »oo N— too 

Proof. Let 

(8-24) flW = At^n^r^Wg, - ^ e n+1 
and 

(8.25) fiW = Av xen M r - iPx N ~ n -^r ln+lX R^ 
so that 

JV-l 

(8.26) R K/ =Y,R {n) - 

n=0 

By the triangle inequality, 

{k^ieiRkA-l^^RkA]) 1 ' 2 

(8.27) 

AT-l 

< (K d E[E[R( n \-L nK )- 1 RW}}) 1 /' 2 . 

n=0 

From the variational formula, 

K'ElR^i-L^r'R^nJ < Av x€nN - n - iPx N - n -^r in+lX £ d n+1 



sup E[RW(-L Qe )~ 1 i?(") 



By the previous lemma, 

(8.28) i d n+1 E[E[R^(-L^ n+ yR^\m en J] < <7M*"V- 
Hence, we have 



oo 

d/2 



(K d E[E[R K/ , (-Ln 1 K )- 1 R K Am t ^ 1 ]]) 1/2 < CMV+W £ ^ 

n=0 

5.29) 

= C(M)r d/2 , 
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with C{M) < oo. Let £ — > oo to complete the proof. □ 



Finally we apply the result to our particular problem. We need to show 
that the difference between (8.1) and 



3.30) 



e d Yj J(Ex)T x E- x v e ^ e (r) kK )W e K dt 

x£Z d /e- 1 Z d ,e 



goes to in P £ -probability as e — > followed by £ — > oo, where K = cq£~ 2 / ( d + 2 ) . 
First of all we would like to replace (8.1) by a term corresponding to the 
first term in (8.22). The difference is 

£ Y J{^x)T x E~ l V eQfi {f)K K ) 

a;eZ<7e- 1 Z< J ) e 

(8-31) 

x [A-VyzM K TyWt - Av xeadmK T x Wlpt >} dt. 
Performing a summation by parts, this can be rewritten as J^Tdt, where 
(8.32) T = e d Y QKA x ) £ ~ lr * W i 

xeZ d /e~ 1 Z d ,e 

and 

QkA x ) = Av yeA' J(£(x + y))u eo>e (T x+ yfjA K ) 



Av yeiZ d nnK p y N) J{e{x + y))u e ^ e {T x+y T]K K ). 



3.33) 

Av yeU d nn K f J y 
Now because V xy fj AK = for x, y e A^-, 

(8.34) (QkAx)t*WH) = r l Av y&TxAl ^ Tx+teKl {Q K/ (x)( Vz - rj y )V yz f). 
Furthermore, 

(QkA x )(Vz ~ VyWyzf) = 2(Qk,£(x)(% ~ %)Vf^yzVf) 

(8.35) 

<2{Q 2 K /x)f) 1 / 2 ((V yz ^f)Y 2 

so 

(Tf) < e^Av^a^-w^Ql^f) 1 ' 2 

(8.36) 

x (r 2 Av yerxA t 7zeTx+ ^(y yz Vf?) 1/2 . 

By the moving particle Lemma 5.2, there is a C < oo such that 

(8.37) r 2 A VyerxAttee A i((y ygy /f) 2 ) < CAv b€TxAeuTx+eeAe ((V b ^J) 2 ), 
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where the average on the right-hand side is over nearest neighbor bonds 
only. Hence, 

(Tf)-s d - 2 D(^f) 

<CAv x ^- l {Q 2 K ^x)f) l ^{Av\ y _ x \ <2 ,(V x ^ x+e ^]f) 1/2 

(8.38) 

< C'Av xeZd/E - lzd>e (Q 2 K)e (x)f) - \e d - 2 D{Jf), 

which vanishes in the limit of small e by the two block estimate (5.6). 

Summarizing the results so far, we have shown that the difference between 
(3.8) and (3.15) vanishes in the limit e J. 0, assuming (8.9). This is proved in 
the next section. Following that, we still have to show that the difference of 
(3.15) and (3.10) is small to obtain the hydrodynamic equation. 

9. Variance estimate for the corrected average current. In this section 
we give the proof of (8.9) which is the key to the renormalization procedure 
of the previous section. 

Lemma 9.1. Let w\ 1 ^\ 2 (m) be as in (8.7). There is a constant C < oo 
such that, for all m G [0, 1] and positive integers K and I with K > 2i, 

(9.1) E[E[{w AlM (m)} 2 \fi n = m]] < Cl~ 2d . 

Proof. Note that (£ X .) A = 1+e l+ ax e° x = e~ A (r/*>A and (( x } x = ^1+°* x 
e~ ax = e A (l - 7] x )\, so that 

(9-2) (6v) Aa = e^ A %, (Ca) Aa = e^ A (1 - fj A ). 

This gives 

(9.3) w AlM - 2m(l - m)(A Aa - A Al ) = B x + B 2 - B 3 , 
where 

(9.4) B X =£aiCa 2 ~ (^a 1 ) Aai (Ca 2 ) Aa2 - Ca^a 2 + (CAi)x Al (U a )x^ 

(9.5) B 2 = (2m(l - m) - r? Al (l - r?A 2 ) - ??a 2 (1 - ?7a 1 ))(Aa 2 - A Al ), 
and, setting cf>(x) = e x — 1 — x, 

(9.6) j B 3 = ?7a 1 (1-% 2 )</'(Aa 2 - A Al ) -j?A 2 (l ~VAi )4>{- C X A 2 -AaJ). 
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Decomposing B\ = B\^\ + + -Bi,3, 

B h 2 = ((e Al >x Al - ^- A ( m ))(a 2 - (a 2 ) AA2 ) 

-((a 2 )A A2 -^- AM )(CA 1 -(CA 1 ) AAl ) 

(9.7) 

+ ((Ca 2 ) Aa2 - (1 - m)e x ^)(^ 1 - <£ Ai ) Aai ) 
" ((Ca,) Aai - (1 - m)e x ^)(U 2 - <a 2 ) AA2 ), 
B 1>3 =me-^[(C A2 - (Ca 2 ) Aa2 ) - (C Al " <C Al >x Al )] 

+ (1 - m)e A H[(f Al - (e" Al ) AAi ) - (Ca 2 - (6v 2 ) AA2 )]. 

We will obtain a bound E [£■[£> 2 |r/oJ] < C£~ 2d for each of Bi tl ,B 1>2 ,B 2 
and B3. Note that .Bi^ is the extra term appearing in 7a 2 — 7 Al . Through 
the proof, C will stand for a finite constant independent of f/Q, though its 
meaning will change from line to line. 

We start with B\^. By the Schwarz inequality, 



E[(£ Al - E[Z Al \fj M ]y(<; Aa - E[( A2 \VA 2 ]y\vn] < Ct 



-2d 



and the same for the analogue of the second term of B\ \. By the equivalence 
of ensembles (see Appendix 2 of [11]), 

-d 



\(&ik Al -E[&i\vM]\<cr*. 

ipute 

E[BlM<Ct 



Hence, it is not hard to compute 

>2 |-i / rm-2d 



When we try to do the same thing for E[B\ 2 \ r ln\i we get terms like 
^[(Ca 2 ~~ -^[Ca 2 |% 2 ]) 4 fe], which are bounded by C£~ 2d , but also a term 
£/[((^Ai) Aa — Tne~^ m ^) \fjo\ for which the same argument does not work. 

Recall (£ai) Aa = ? ?AiC~ Aai and rewrite the difference as two terms 

[??A 1 e- A( ^ A i ) - me~ A(m) ] + [^-^^-(Aa^a^-a^)) _ ^ 

The function pe~ x ^ is bounded and Lipschitz and E[E[(t/Ai — w) 4 |^n = 
m]] < C£~ 2d . Hence, for the first term, we have the required bound 
E[E[(??A 1 e- A ^ A i ) - me- A ( m )) 4 |%} =m\] < C£~ 2d . For the second term, note 
that ?7aiC~ a ^ Ai ^ is uniformly bounded. For |x| < A, \e x — 1| < C\x\ so 
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E^IOAa^aJ-A^aJI <^)(e~^ A i^ A i)- A ^ A i))-l) 4 |r/ Q ] is bounded above 
by a constant multiple of E[£'[(Aa ! (%J — A(% i )) 4 |r/n]]. By Lemma 10.3, 

(9-8) E[(A Al (%J-A(^)) 4 ]<Cr M 

for i = 1,2. On the other hand, r/Ai e ~^ AAl ^ Al ^ and ??Ai e ~ A ^ Al ^ are both uni- 
formly bounded and by Chebyshev's inequality and (9.8), E[i£[:[L(|A Al (r7 Al ) — 
A(t? Ai )| > A)\fjn] < CA- 4 l~ 2d . We conclude that 

E[E[Bl 2 \fj n ]]<Ce- 2d . 

Turning to B2, it is not hard to see that 

E[£[((2m(l - m) - t? Ai (1 - m 2 ) - fj Aa (l - ?? Al ))) 4 fc]] < C£~ 2d . 

We claim that 

(9.9) E[E[(\A 2 -\A 1 f\vn]}<Ce- 2d . 

Then by Schwarz's inequality, we conclude that EfE'f.Bf 1%}]] < C£~ 2d . To 
prove (9.9), rewrite A Aa - A Al = Aa 2 (?7a 2 ) - A Al (r/ Al ) as 

[A Aa (?? A2 ) - A(r? A2 )] + [A(t7a 2 ) - A(t? Ai )] - [A Al (r/ Al ) - A(t? Ai )]. 

The first and third terms can be handled by (9.8). Using the fact that the 
variance of r\ x is bounded by a constant times 1/A', it is not hard to check 
that also E[£[(A(% 2 ) - A(%J) 4 fe]] < Ci~ 2d , proving (9.9). 

Finally we consider B3. Let X = A Aa — A Al . From (9.9), we know that 
E[E[X 4 \fj n ]] < C£~ 2d , so it will suffice to bound E[E[B$ \fja\] < CE[E[X 4 \fj a ]]. 
There is a constant C = C(A) < 00 so that if \x\ < A, then |</>(x)| < Cx 2 . So 
1{\X\ < A)cf> 2 (X) < CX 4 and, hence, 

E[£[l(|X| < A)B 2 \f)n]] < CE[£[X 4 fe]]. 

On the other hand, it is not hard to check that 

y = r/ Al (l-% 2 )(e AA 2- AA i -l)-77 A2 (l-7/ Al )(e AA i- AA 2 - 1) 

is uniformly bounded. Hence, 

E[E[t{\X\ > A)Y 2 \fj n ]] < CE[E[t(\X\ > A)\fj n ]] 

and by Chebyshev's inequality, the last term is bounded above by A _4 E[i?LY 4 | 
fjn]]. Finally, E[E[1(\X\ >A)X 2 \fjci}} <2E[E[1(\X\ > ^)fe]] + 2E[£[X 4 |r? n ]]. 
The first part is bounded by 2A~ 4 E[£'[X 4 |?/n]] by Chebyshev's inequality 
again. □ 

Lemma 9.2. 

K£ai>a a -E[hAfj Al ]\<cr d . 
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l+e L 



Proof. Fix A and let P denote the product measure with P(rj x = 1) 
and E the corresponding expectation. We claim first that 



(9.10) 
Now 

(Cs 



|A| 



P(12x€AVx = M) |A|-M + 1 



e~ x E 



Ca| 

xeA 



M 



E 



Ca| J2 r i*= M 



xeA 



3 Aa 



^1 1 + e°«+AA 
P A\ 

E — — - 

^1 1 + e«-+AA 



-P(E, gA - W % = M-l) 



1 - P(r/ X . = 1) + Pfe = 0) 



P{E v eA-{x}Vy = M-l) 



Now 



(9.11) 



P(E y eA-{x}Vy = M-l) 



+ 



P(E v eA-{*}Vy = M) 



By (9.10), 



(9.12) 



1 + e «,+A A 1 + e a x +X A P{E y& A-{x} Vy = M - 1)' 



P(£ v eA-{x}Vv = M) 



P(ZyeA- {x} Vy = M-l) 
f IAI-1 



[ |A| - M 
By the inductive hypothesis, 

E 



e- XA E 



Ca-{x}\ E % = M 
yeA-{x} 



(9.13) 



<A-{x}| E % = m 

yGA-{x} 



(Ca-m>a a + 



c 



{*} |A|-1 

-w l A l- M i C 
|A|-1 + |A|-1 



□ 
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10. Fick's law. Our goal in this section is to prove that, for any smooth 

[ T e d J(ex)[e- X T X & K - D^M 1 ^)^)' 1 

x(£e- 1 Z d /Z d 

(10.1) 

x W\ £ /-i -ff 1 ** 1 i )}dt 

v 'X+02S 1 e <x—02E i e /J 

vanishes in P e probability, as 81,62 — > 0. Q e K is given is (3.28). 

Let mi and r?72 be the particle densities on = Ak and K 2 K = tkc^k- 
Let Aj(m) = X\(m), i = 1,2, be the empirical chemical potential on A l K as 
in (8.4). Let 

The expectation is with respect to the product measure on Ai with chemical 
potential Ai(mi) and on A2 with chemical potential A2: 

(10.3) Z~ x exp< ^2{a x + Xi(mi))r] x + ^2(a y + X2(m2))r]y>. 

IxeAi yeA 2 ) 

The function F(mi,rri2), which depends on the densities mi and 7772, as well 
as on the field configurations on the blocks A]^ and A^, is given explicitly 
by 

F(m u m 2 ) = /^[(^("^-Mmi) _ i) mi (i_ m2 ) 

(10-4) 

- ( e M™i)-A2(m 2 ) _ l) m2 (l _ mi )]. 

Consider as well the variant of F where we use the annealed chemical po- 
tential instead of the empirical chemical potential 

F(mi, 777 2 ) = if- 1 [( e A (" l 2)-A(mi) _ 1)mi ^ _ m2) 

(10.5) 

- ( e A(mi)-A(m 2 ) _ i} m2 (i _ mi )]. 

For mi 7^ m,2, consider the quotient $(7771,7772) = F(mi, m2)/{K~ l {m2 — 
mi)). From the boundedness of the field a, it is bounded above and be- 
low away from and Lipschitz. When mi = 7772 = m, we have $(777,777) = 
2A'(?77)777(1 — 777). We call (7(7771,7772) = 1/$ (777-1,7772). G defined in this way 
is also bounded and uniformly Lipschitz on [0, 1] x [0, 1]. 

We will prove (10.1) is several steps. The first thing we want to do is 
replace the term 

(10.6) Av x J(ex)e- 1 T x Q e K = Av x e~ x J{ex)T x v(f) kK )W e K 
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in (10.1) by 

(10.7) Av x e- X T x TWfc, V = [Av y p 5 J' 5 * J(ey)r yJ D eoe (^ 1 ^ 1 )]G, 

where the weights p y 1,&2 are defined on the convex hull of two boxes of 
side length e^di whose centers are separated by the vector e~ 1 52e and 
J^yPy 1 ^ 2 = 1- Here v(m) = a eoe (m)/m(l — m), G = G(mi,m 2 ) and we have 
abused our definitions mildly by writing r x J{ey) = J(e(x + y)). 

The weights p y 1,&2 are defined as follows. We can assume without loss 
of generality that K divides both e~ x 8\ and e~ 1 5 2 evenly. Divide the box 
centered at the origin of side length e~ 1 5\ into boxes of side length K and 
label their centers f3. For any a y , y £ 1> d , let 

(10.8) Av y p s y ^a y = Avp(2S 2 )- 1 sK £ a PHl _ l/2)Ke . 

i=-S 2 e- 1 K- 1 

This average has the property that if a y — Wy+Ke/2 — ^y—Ke/2 , > 

{eK)^Av yP ^a y = (25 2 )-\f 1 l\ e £ -l e ~ V%\ e ). 
A summation by parts gives 

£ d-i K -iJ2r x TF = e^Ji^TxDeUVo 16 ' 1 ) 

X X 

which is the right-hand side of (10.1). In other words, if we can replace v by 
T as described in (10.6) to (10.7), and then replace Wk by K~ l F, we end 
up with the right-hand side of (10.1). 

We start with (10.6) to (10.7). The definition of constants C will change 
from line to line, but will always denote a finite constant independent of the 
parameters e, 5\ and 5 2 . Let \l/ = Jv{rjQ K ) — T. We have 

(10.9) (W K *f) = K- 1 Av x ^ yeK% {(V xy r 1 ){V xy yf)). 

Now V xy ^ f = X y f + f (T xy r])V xy ^> . The first piece ^V xy f is rather stan- 
dard. Write V xy f = Vt y ^fV xy Vf, where V+ / = / (T^rj) + / ( V ) . Changing 
variables T xy rj^rj, we can write {W xy r}V xy f) = ((V+ / ^) v / /V X j / r/V xy s/J). 
By Schwarz's inequality and the moving particle lemma, for any q > 0, 

Av xeA i Aj > k e*V xy r,V xy f) < Cq(* 2 f) 

(10.10) 

+ K 2 q~ 1 Yl ((V W V7) 2 ). 
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In the second piece, f \T xy rj)V xy ^> , 

V^* = GiT^V^Avup^ JD eoe (77* ie_1 ) 

(10.11) 

+ Av u p%> Sa JD eoe (ffc e ~ )V xy G. 

From the definition (10.8) of the weig hts in the first term above, there 

are only a finite number of terms in the sum for which V xy D eQe {f)^ £ ) is 
nonzero and those that do give a term of the form GJ(62)~ 1 £K(D eoe (fjf l 1£ 1 ± 
(6\)^ 1 e) - D eoe (fj^ 1£ For fixed 61,62 > 0, this is o(eK). As for the second 
part, since G is uniformly Lipschitz, |V xy G| < CK~ d . So we have 

(10.12) Av xeAlKtyeA 2 K (fV xyr] V xy *) = o(eK). 
Hence, taking q = eK in (10.10), 

£~ X Av x&d/e -i zd (T x mW K f) - e- 2 Av Xje ((V xx+e ^) 2 ) 

< CAv x ((T x ^) 2 f) - \e- 2 Av Xie {{V xx+ eVf?)+o{l). 

The first terms on the right-hand side vanish as e — ► followed by 6\ — > 
and 62 — ► by the two block estimate. To conclude, we have 

(10.13) limsupsupfe"^^*^/) - e~ 2 Av x , e ((V xx+e ^f) 2 )} = 0, 

$2,Sl,S / 

which shows that we can replace (10.6) by (10.7). 

The next step is to replace Wk by F{pi\, 771,2). We do it in two steps: First 
replace Wk by F (Lemma 10.2) and then F by F (Lemma 10.4). Before 
proving Lemma 10.2, we need a preliminary estimate which shows that the 
integration by parts property of our long jump currents is preserved under 
conditional expectations. 

Lemma 10.1. Let Q = AU B be subsets ofZ d with max x ^A,yeB \y — x\ < 
K . Let x £ A and y G B and w xy be the current given by (3.1). Then there 
is a C < 00 so that 

(E[w xy \fj A ]f) < CK ]T ((V M ,V/) 2 ), 
ueA 

(E[w xy \fj B ]f) <CKj2((V xv s/f) 2 ), 
{E[w xy \fj A ,fj B ]f)<CK Yl ((V^v/J) 2 )- 

u£A,v£B 

Proof. Integrating by parts, 

(10.14) (E[w xy \fj A ]f) = (( Vy - Vx )V xy E[f\fj A ]). 
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There exist bounded a u , u£ A so that 

(10.15) VxyE[f\fj A ] = E E[a u V uy f\fj A }. 

ueA 

In fact, 

P(v u \E z&A Vz + l) , n m s P(v u \E z& aVz-1) 



a, 



Vyi 1 ~ Vu) K ' T' £A " ' + (1 " Vy)Vu- 



where rj u is the configuration changed only at u. Hence, 

(10.16) {E[w xy \fj A ]f) = J2 (KVuyf), 

u&A 

where b u = a u E[(r]y — r] x )\f]A\- The result for i^w^l^] follows by Schwarz's 
inequality. The other results are proved in exactly the same way. □ 

Lemma 10.2. Let 

(10.17) n 1 =Av x T x T[W K -F(m l ,m 2 )}. 
For any 7 > 0, 

lim sup lim sup lim sup sup{e _1 (/fii ) 

(10.18) 

- je- 2 Av x>e ((V X:X+e ^f) 2 )} < 0. 

Proof. Let F 1 = E[W K \mi] and F 2 = E[W K \m 2 ]. By Lemma 10.1, W K , 
Fx, F-2 and F all satisfy (3.22). From these estimates and the two block 
estimate, we can replace V by T, where T(j)(ey)D eoe (rj K ))G(mi,m 2 ). We will 
write fii as K' 1 Av x Tt x [(W k -Fi~F 2 + F) + (F l + F 2 - 2F)\. Now 

W K - F x - F 2 + F = (£ Al - E[£ Al \fj Al })(( A2 - E[( A2 \VA 2 ]) 

- (CAj - ^[CAiI%i])(^A 2 - E [£A 2 \fjA 2 ])- 

Hence, 

((W K -Fx-F 2 + F) 2 ) = 0(K- 2d ). 
Now Fi + F 2 — 2F consist of a term of the form 

(10.19) (f Al - E[i Al \fj Al ])E[CA 2 \fjA 2 ] - (6v 2 - E[Z A2 \fj A2 ])E[( Al \fj Al ] 

and another one with the roles of £ and 77 reversed. Summing by parts, we 
can rewrite Av x Tr x [Fx + F 2 — 2.F] as 

(10.20) A^esXTXe - T- Ke )j(U - E[£ A |r? A ]), 
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where j = /i(t7a).E[Ca|^a]> plus an analogous term with £ replaced by £. Since 
we have 

(10.21) E[((r Ke - r. Ke )j) 2 (U - E[U\fj A }) 2 } = 0(K~ 2d ), 
we have the lemma. □ 

Now we finish the proof with the replacement of F by F (Lemma 10.4). 
Note here the special role played by the scale K = 0{e 2 l d+2 ). We need first 
a preliminary result on the difference of the chemical potential and the 
empirical chemical potential. 

Lemma 10.3. Let a x , x^7L d be a random field taking values in [—B,B] 
for some B < oo and satisfying the mixing conditions (2.1). Let Xx(m) be 
the empirical chemical potential on a block Ak of side length K , given by 
(8.4). Let A(m) be the annealed chemical potential given by (2.7). Let 7 be 
as in (2.1). There is a constant C independent of m G (0, 1) so that 

(10.22) E[\X K (m) - A(m)H < CK^ d l 2 . 

Proof. Let L m denote the real interval [log — B, log + B] . Since 

\a\ < B, we have A(m), Xk(tu) G I m . Recall the definition (8.4) of the empir- 
ical chemical potential and denote p(x) = e x /(l + e x ). By Taylor's theorem, 

m = Av xe A K p(a x + A) + (X K - X)Av x(LkK p(a x + A)(l - p{a x + A)), 

for some A € I m . Therefore, 

E[|A* - A|^] < r^E[\Av X £A K {p{oc x + A) - m)\\ 

where r m = inf 3/g / m p 7 (y)(l — p(y)) 7 . By the mixing conditions, 

E[\Av xeAK (p(a x + A) - m)| 7 ] < CK-^ 2 E[\p(a + A) - m| 7 ]. 

By Taylor's theorem, 

p(ao + A) — m = (ao — a*)p{a + A)(l — p(a + A)), 

for some |a| < B, where \a*\ < B is chosen so that p(a* — A) = m. Therefore, 

E[\p(a + X)-m\y]<4B^R m , 

where R m = sup yeIm p J (y)(l -p(y)) 7 . One can check sup m6 ( 01) R m /r m < 00, 
and, hence, (10.22) follows. □ 
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Lemma 10.4. Let 
(10.23) n 2 = Av x T x T[F(m 1 ,m 2 ) - F(m 1 ,m 2 )}. 

Then for any 7 > 0, 

lim sup lim sup lim sup sup{e _1 (f£l 2 ) 

(10.24) 

- 7e- 2 Av X)e {(y x , x+e ^f) 2 )} < 0. 

PROOF. Write F-F = K' 1 [A ± + A 2 ] with 
(10.25) A 1 = l(\ 2 - A)(m 2 ) - (A x - X)(mi)]H, 

where H = H(mi,m 2 ) = [m ie -^ mi \l-m 2 )e x ^ +m 2 e- x ( m2 \l-m 1 )e x ^} 
and 

A 2 = (t>(U)b(mi,m 2 ) - 4>(-U)b(m 2 ,mi), 

(10.26) 

U = (X 2 -X)(m 2 )-(X 1 -X)(m 1 ), 

where b(x,y) = xe~ x ( x '(l — y)e x ^ y > and 4>(x) = e x — 1 — x. Note that |A — A| < 
2B always and, hence, we have 4>(x), <p(— x) < Cx 2 for relevant x. It is also 
easy to check that b is uniformly bounded. Hence, we can estimate 

(10.27) e- 1 K~ 1 Av x T x TA 2 < Ce^K^Av^Xr - A| 2 , 

where Xk = Xk{ti K ) and A = X(fj K ). Next we turn to the A\ term. Summing 
by parts, 

(10.28) e^K^Av^TH^TKe - I)(X K - A)] = a' 1 K~ x Av x t x B x (X k ~ X). 

B = (r-Ke - I)TH. Let a > 1 be as in (2.1) and a* = a/(a — 1) be the 
conjugate exponent. The right-hand side is bounded above by 

(10.29) Cq a *Av x T x \Bf + q- a (eK)- a Av x T x \X K - X\ a . 
By the two block estimate, for any C > 0, 

(10.30) limsu V sup{CAv x (T x \Bff}-e- 2 Av x , e ((V xx+e ^f) 2 )}<0. 

6^0 f 

By Lemma 10.3, and since K > co£~ 2 ^ d+2 , there exists Co < 00 such that 

(10.31) limsupEfKeK)" 1 ^ - A)| 7 ] < C . 

By the mixing conditions, 

n{Av x r x {(eKy l (X - X)) a - E[((eK)~ 1 (A - A)) a ]) 7/a ] 

(10.32) 

< (7 e 7tZ 2 /2a(d+2) 
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which is summable in e" 1 = L = 1, 2, . . . as long as 7 > 2(d + 2) jd 2 and a 
is sufficiently close to 1. By Chebyshev's inequality and the Borel-Cantelli 
lemma, for almost every realization of the random field, 

(10.33) limsup q- a Av x T x {{eK)- l {\ - \){m K )) a < q- a C . 

Letting q — > 00 completes the proof for A±. The same argument with 7 = 2 
and the estimate (10.27) show that e~ x K Av x t x ^ Ai — ► with probability 
one. □ 

11. Continuity of the diffusion coefficient. We introduce a notion of reg- 
ularity on D. There is a finite C such that, for any x and y in [0, 1], 

(11.1) \D(x) - D{y)\ 2 < C(x(l - x))' 1 ^ - y\. 



Lemma 11.1. The diffusion coefficient D(m) satisfies (11.1) 
PROOF. Fix a vector (3 G R d and let 



(11.2) F(m) = \l m{l - m)f3 ■ D- l {m)(5. 



From the proof of the hydrodynamic limit, we have the following represen- 
tation of the diffusion coefficient: 

(11.3) F 2 {m) = C lim K~ d E[V K ], 

K^oo 

where 

(11.4) V K =supj2]T/3- e ((Vy-r] x )V xy f) m -V K . m (f)\, 

y€T Ke A K 

where T>K.m is the Dirichlet form on Ak Urxe^-K- The expectation is with re- 
spect to the ergodic invariant measure on that box with density m. The con- 
stant C comes from the nonstandard average as described in the Introduction. 
The exact value of C is not relevant. 

In order to compare two densities m and m + h, we will produce a coupled 
measure. Independently, at each site x, place a red particle with probability 

1+eax +x( m ) > a g reen particle with probability 1+eCtx+X{m+h) - 1+eax+ x [m) and 
no particle with probability 1+f , ax l\( m+ h) ■ Let p x 6 {0, 1} denote the presence 
or absence of a red particle and ^ x the same for a green particle. p x + r y x = r} x 
is the whole configuration which has measure (2.6) with average density 
m + h. We can also couple the dynamics as follows. The red particles have 
priority in the sense that when a red particle tries to jump on top of green 
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particle, the two particles switch positions. Otherwise, the particles evolve 
as before. In the common usage, red are first class particles and green are 
second class particles. If we are colorblind and see only the total particles, 
the evolution of r](t) is as usual. If we cannot distinguish green particles 
from empty sites, the evolution of p(t) is also as usual. Hence, if Wf^ = 
&^K,y&T Ke A K w§. y , where 

p x (l - p y )(l + e a *~ a y) - p y (l - p x )(l + e°y- a *) 
is the red current, we can write 

(11.5) F\m) = \\Wq 2 m>m+h = lim K- d E[V p k ], 



W 



xy 



where 
(11.6) 



V£ = SUp{2(/? • W^f) m ,m+h ~ ^K,m.m+h(f)}, 
f 



where T>K,m,m+h is the Dirichlet form for the coupled process. Similarly, 
(11.3) and (11.4) can be rewritten in this language as 



(11.7) 

where 
(11.8) 



F 2 (m + h) = \\W\\ 2 m , m+h = lim K~ d E[V k ] 

K— >oo 



V fc = SUp{2(/3 • W K f)m,m + h ~ ?>K (/)}• 
/ 



Let Vfn = W-W.We have 
(11.9) \F(m + h) -F{m) 



\W\ 



\W P \\\ < \\W~t\\ 



and 
(11.10) 

where 
(11.11) 



\W p \\ 2 m , m+h = lim K~ d E[Vl], 



VJ = SU P {2(/3 ■ Wlf) m ,m+h ~ V K ,m,m+h (/)}• 
/ 



Now 



(P ■ Wlf) 



KJ I ra,m+h 



(P ■ W K f) m ,m+h ~ (P ■ W£f) 
= K' 1 J2P' e Av xeA K ((ly ~ lxWxyf)m,m+h- 

e y&TKe^K 

By Schwarz's inequality, 

((ly-lx)V xy f) < ((ly-lx?) 1/2 ((Vxyf) 2 ) 1/2 

<Ch^ 2 ((V xy f) 2 ) 1 / 2 . 



(11.12) 
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From the moving particles lemma (5.2), we conclude that 

(11.13) \F{m + h) -F(m)\ <CVh. 
From the definition (11.2) of F, this gives 

y/m(l -m)\yfp- D~ l (m + h)(3-^J(3- D~ 1 (m)/3| 

< Ch 1 / 2 + |/3 • D- l {m)f3\ 1 l 2 \{{m + h)(l - (m + h))) 1/2 - (m(l - m)) 1/2 |. 
Since 5" 1 1 < D < 51 for some finite this implies that 

(11.14) \D(m + h)- D(m)\ < C(m(l - m))" 1 / 2 ^ 1 / 2 . 

By particle-hole duality (77 1 — >• 1 — 77 and a t— > —a), we conclude that \D(m — 
h) - D(m)\ < C(m(l - m))- l l 2 h 1 / 2 as well. □ 

Equation (11.1) say that D{m) is Holder 1/2 in (0, 1), but it says nothing 
about continuity at the end points. We prove it separately now. 

Theorem 5. The diffusion coefficient D(m) is continuous on [0,1]. 

Proof. By the symmetry (r/,a) 1 — >- (1 — 77, —a), it follows that D{m) = 
D(l — m), where D is the diffusion coefficient of the dynamics with a re- 
placed by —a. Hence, one only needs to check the continuity at 0. 

For any subset A S Z rf , let 



(11.15) 



VA 



rife 



where p x 



e a x +X(m) 
I _|_ e a x +X(m) ' 



These form an unnormalized orthogonal basis of L 2 (u m ), where u m is the 
infinite product measure (2.6) with density m. Let H n denote the subspace 
spanned by fjA with \A\ =n. Let 



(11.16) 



Q2 = (J H n- 



n>2 



We can write the variational formula (2.12) as 



(11.17) (P, a(m)(3) = 2 inf < E 

U(a) 



H 6 e((% -Vo) 



r(m) 



where 
(11.18) 



P e + r_ e U - U, 



(11.19) r(m) = supE 



26 e ( (r]e - Vo) Voe r x g \ - ( V 0e Y T *9 
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Setting g = shows T(m) > 0. Hence, 



(11.20) (j3, a(m)(3) < 2 inf E 

U(a) 



J2((Ve-V0) 2 m + T- e U-U)' 



We now work toward a bound in the other direction. For any /, 

(11-21) (fae-»fc)Voe/) = <uW) 

and WQ te € Hq U Hi U i?2 with projection onto H2 of the form cfjo >e . Let 
5 = J2\a\>2 (JAVA- Note that since g is local, this is a finite sum. Furthermore, 
note that T(m) only acts on J2 X T x9i where the sum is over a box large enough 
that Vo te T x g = for any x in its exterior. Because of the sum over shifts, 
we can assume without loss of generality that g = go,eVo,e + /> where f £ G2 
and / = E\A\>2 A ^ x+e} for any JAVA- Then 



(11.22) ( (Ve-Vo)^0,eJ2 Tx 9) = -90,e(Pe - Po){(Ve - r]Q 



and 



(11.23) V 0e ^r^ = ^ + Vo, e ^r x /, 

X X 

ie = {Ve ~ Vo)(fj-eTe90,e ~ (Pe ~ Vo)h,e + V2eT-e90,e 



(11.24) 

Hence, 

(11.25) r(m) = supi -2E 



90, e 



^2b e ((Ve-r)o) 2 )(pi -Po)9\ 



Oc 



Vl'(m) 



(11.26) ¥(ro)=iafE 



E( Ke+ V 0,eE T -/ 



The supremum is over G2 = {/ 6 G2 '■ f -L span^a+e, x G 

We claim that there is a Ci > such that, for sufficiently small m > 0, 

(11.27) ^(m)>Cim 2 E[g2 e ]. 
To prove (11.27), note first that 

(11.28) Z e = il + i\ + (r? e - m ) {Pe - Po)SO,e, 

where E -f/j. G2 is preserved by shifts, and for f £ G2, ((rj e — i]o)^o,e J2x T x x 
/) = 0. Hence, 



*(m) =E 



(11.29) 



E^^-^o) 2 )^-^) 2 ^ 

E^^+d-Vo.eE^/ 



+ inf E 

/6G 2 
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NOW £ 2 = (fj e - fj ){fj- e T e go,e + ??2eT_ e 50,e) and £ = (Pe ~ Po) (?7-eT e 30, e + 

fi2e T -e9o,e) ■ One can check then that there is a c\ < oo such that 



(11.30) E 



X)((??e-??o) 2 )(Pe-Po) 2 3o, f 



+ E 



< 



Cl m 3 E[5 2 e ]. 



Hence, it is enough to show that, for some fixed e, there is a C2 > such 
that, for small enough m > 0, 



(11.31) 



inf E 

/€G 2 



V. x J 



2\ n 



>C 2 m 2 E[£ 2 ]. 



Without loss of generality, we can write / = X)a/a^7a, where G 4 and 
^ r neJ 4 for n < 0. Now one can check using the explicit form of £ 2 that 
Vo,eT~xVA -L f° r an such ^4 except A G {{0, 2e}, {0, e,2e}}. So the infimum 
in (11.31) is achieved with / = /o,2e%,2e + /o,e,2e?7o,e,2e- Then Vo^E^J is 
given by 



(11.32) 



(/0,2e??2e ~ T_ e /o,2e??3e ~ T e fo,2eV-e + 7"2e/o,2e??-2e 
- (Pe ~ Po)fo,e,2eV2e ~ {Pe ~ Po)T e fo,e,2eV-e 

+ T~2efo,e,2eV-2e,-e + T"-e/o,e,2e%e,3e) 
X (??e-??0+Pe-Po) 

and we can compute E[((£ 2 + — Vo, e r a;/) 2 )] explicitly to get 
E[{p e +p - 2p e po} 

X {P2e(l - P2e){-r- e g0,e ~ h,2e ~ (Pe ~ Po)h,e,2ef 

+ P-e(l -P-e)(T e go,e ~ T e fo,2e ~ (Pe ~ PoK jo,e,2e) 2 
+ p 3e (l -p 3e )(r_ e / ,2e) 2 +P-2e(l -P-2e)(T2e/o,2 e ) 2 
+ p_2e(l ~ P-2e)P-e(l ~ P-e)(T2efo,e,2e) 2 

+ P2e(l -P2e)P3e(l -P3e)(T- e /o, ei 2 e ) 2 }]- 

Using the bound on the field a, one has p x > C^m for some C3 > 0. Then it is 
not hard to check that this quadratic form is bounded below by C47b 2 E[<7q J 
for some C4 > 0. Using the upper bound (11.30) on E[V^ e ((£g) 2 )], we obtain 
(11.27). 

Now (11.27) implies that 



[11.33) 



r(m) < C 2 m 2 E 



HZ 
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Recalling (11.17), we have therefore shown that 



(11.34) (P, a(m)(3) > 2 inf E 

U(a) 



EK(^ " Vof) ~ C 2 m 2 ](f3 e + T- e U - Uf 



Together with the lower bound (11.20) and the formula (2.13) for -D(O), we 
only have to prove that, for any C bounded, 



(11.35) 



lim A'(m)2 inf E 

m-»0 U(a) 



2z~ l inf E 



Et((^ " %) 2 ) + Cm 2 ](J3 e + r. e U - Uf 



J2(e ao +e ae )(Pe + T- e U-Uf 



Using ((rj e — r?o) 2 )m = Pe + Po ~ %PePo and the form (11.15) of p x , we have 

e A(m)( e ao +e a e ) 



(11.36) 



<(»7e - J/o)' 



(1 _)_ e a o+^( m ))(l + e Qe +^( m )) ' 



Since — B < on < B, 

e A(m)( e ao +e ae) 



(11.37) 



(1 + e -B+A(m))(l + e B+A ( m )) 



< (foe - Vof 



=A(m) f„ao 



< 



(e a ° + e Qe ) 



(1 + e -B+A(m))(i + e -B+\(m) ) " 

It is elementary to check that as m — > 0, A(m) ~ logm — log z, and A'(m) ~ 
m. Hence, (11.35) follows. □ 

12. Uniqueness. We have the following theorem of uniqueness of weak 
solutions. 

Theorem 6. Let u\ and u 2 be two weak solutions of dtu = V • D{u)Vu 
[see (2.15)] satisfying (2.14), with the same initial data and suppose that D 
satisfies (11.1). Then u\=u 2 . 

Proof. Taking in (2.15), 4> e to be the solution at time e > of the 
standard heat equation on T d with initial data Sq, and denoting by f £ the 
convolution of a function / with (j) E , we have, for i = 1, 2, 

(12.1) di(0 £ = V-(£>(OVui) £ 

as an equality of smooth functions. Let ips( x ) be an approximation of \x\ 
with ipf(x) = (2vr5)- 1 / 2 exp{2; 2 /25}. From (12.1), we have 

(■T 



JT d 



lj)'s({ui - U 2 ) £ ){V(u 1 - U 2 )) e 

x (D(ui)Vui - D(u 2 )Vu 2 ) £ d8dt. 
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Let s — ► 0. Since Ui are bounded, F((ui — u 2 ) £ ) — ► F(u\ — u 2 ) in L°°. Prom (2.14), 
Vui a 2 = L 2 ([0, T] x T d ) and, since D is bounded, D( Ui )Vui £ L 2 as well. 
Hence, (Vuj) e — > Viij and (Z?(nj)V?Xi) e — > D(ui)Vui in L 2 and we can take 
the limit of the above formula to obtain 

! M(ui-u 2 )(T))d6 
Jj d 

= ~ I I -u 2 )(Vui - Vu 2 )(D( Ul )Vui - D(u 2 )Vu 2 )d9dt 

JO JT d 

= ~ I \ i>&{ui-u 2 )(Vui-Vu 2 )D{ui){Vui-Vu 2 )d0dt 
Jo Ji d 

+ f T I ips{u l -U2){Vu l -Vu 2 ){D(u 1 )- D(u 2 ))Vu 2 d6dt. 
Jo Jf d 

By Schwarz's inequality and the fact that cl < D, we obtain from the pre- 
vious line that 

/ M(ui-u 2 )(T))d8<C [ T [ ^(tti-t i2 )|D(u 1 )-i)(« 2 )| 2 |VM2| 2 J<it 
Jj d Jo Jj d 

for some finite C. From (11.1), this is bounded above by 

(12.2) C [ T [ W>( Ul - U2 )\ Ul -u 2 \)( ^ U2? )d0dt, 

Jo JT d \u 2 (l-u 2 )J 

with perhaps a new finite C. Let F = J^i2u 2 ) £ L 1 by (2.14). Let v = u 1 — 
u 2 and h(x) = (2-7r)~ 1 / 2 |x| exp{— x 2 /2}. Note that h is uniformly bounded 
with h(0) =0. Rewrite (12.2) as C J T f Td h(6- 1 / 2 v)F dO dt. The integrand is 
dominated by H^Hoo-? 1 £ L 1 and h(5~ 1 / 2 v ) — > pointwise. Letting 5 — > by 
the dominated convergence theorem, we obtain 

(12.3) hm/ T / (^(m-n 2 )|m-n 2 |)f )d9dt = 0. 
s^oJo Jr d \u 2 {l-u 2 )J 

One the other hand, by the monotone convergence theorem, we have 

(12.4) lim / ^ s ((u 1 -u 2 )(T))d0= [ \m - u 2 \(T) d9, 

S^0jf d Jjd 

from which we conclude that Jj d \u± — u 2 \{T) d9 < 0. Hence, ui = u 2 . □ 
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